Fermi-Edge Singularities in the Mesoscopic X-Ray Edge Problem 
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We study the x-ray edge problem for a chaotic quantum dot or nanoparticle displaying mesoscopic 
fluctuations. In the bulk, x-ray physics is known to produce Fermi edge singularities - deviations 
from the naively expected photoabsorption cross section in the form of a peaked or rounded edge. 
For a coherent system with chaotic dynamics, we find substantial changes; in particular, a photoab- 
sorption cross section showing a rounded edge in the bulk will change to a slightly peaked edge on 
average as the system size is reduced to a mesoscopic (coherent) scale. 
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The x-ray edge problem for metals refers to a singular- 
ity at threshold in the photoabsorption spectra associated 
with the excitation of a core electron to the conduction 
band by an x ray. The singularity takes the form of 
a peaked or rounded edege and was studied intensively 
over the last four decades using different methods 
The sustained attention devoted to this problem devel- 
oped because the basic phyiscs originally studied for x- 
ray absorption recurs in many other contexts. Recently, 
photoabsorption in semiconductor quantum wellsj^] and 
electrical conduction through localized states Q have 
been particularly active areas of Fermi edge research. 
Motivated by the growing interest in nanoscale systems 
like quantum dots or nanoparticles H, |M 111 i we address 
Fermi edge physics for such confined geometries, assum- 
ing generic chaotic dynamics of the electrons. 

For bulk metals or degenerately doped semiconductors, 
the photoabsorption cross section A{uj) near the thresh- 
old energy Wth involves two competing effects, Mahan's 
enhancement and Anderson's orthogonality catastrophe 
(AOC). Its form is Q 

Here 6i is the partial-wave phase shift at the Fermi en- 
ergy associated with the core hole potential Vc for orbital 
channel /; lo is the optically excited channel. The phase 
shifts obey the Friedel sum rule Z — X); 2(2/ -I- 1)Si/tt 
with Z = —1 here. AOC |3| reduces the overlap between 
ground states before and after Vc is applied, leading to a 
rounded edge. On the other hand, the sudden presence 
of a charge associated with the missing core electron en- 
tails a many body response that enhances the threshold 
cross section, resulting in a peaked edge. Whereas all I 
channels are sensitive to Vc and so contribute to AOC 
(second term in the exponent), only the optically excited 
lo experiences the many-body interaction (first term). 

The many-body enhancement depends, via the dipole 
selection rule, on the symmetry of both the core elec- 
tron and local conduction electron wavefunctions. We 
will assume the latter to be of s-type, and distinguish 



between core electrons with s-symmetry (K-shell) and p- 
symmetry (L2,3-shell) 0, referring to the photoabsorp- 
tion threshold as a K- or L-edge, respectively. Further- 
more, we take the core hole potential to be spherically 
symmetric: Si = for 1^0. The Friedel sum rule 
then implies Sq = — 7r/2, so that either a rounded K- 
edge {lo = 1, photoabsorption dominated by AOC) or 
a peaked L-edge {lo — 0) is found. For mesoscopic sys- 
tems, in addition to So ~ — 7r/2, we shall consider smaller 
values for situations where other charges (e.g. on nearby 
gates in a semiconductor heterostructure) help to screen 
the core hole. 

Compared to bulk metals, some characteristics intrin- 
sic to nanoscale objects will significantly alter the x-ray 
edge properties. First, the fact that both the number of 
electrons and the total number of levels are finite will 
lead to an incomplete AOC and many-body enhance- 
ment. Furthermore, the lack of rotational symmetry 
eliminates independent orbital momentum channels. Fi- 
nally, mesoscopic fiuctuations will affect both energy lev- 
els and wavefunctions. In tackling these issues, the key 
ingredient is the joint statistics of the electronic proper- 
ties before and after the core hole excitation j^. 

The influence of these effects on the average and fluctu- 
ations of the Anderson overlap was modeled using para- 
metric random matrices and for disordered systems 
jllj . These studies mainly addressed the regime where 
the added potential can be treated perturbatively, as, for 
instance, for scrambling phenomena in Coulomb blockade 
experiments. Here, we study first AOC and then absorp- 
tion spectra near threshold for a model corresponding to 
the x-ray edge problem. In this case the perturbing po- 
tential Vc associated with the core hole has a very short 
range, but on the other hand strong perturbations are 
physically relevant. Possible experimental realizations 
are briefly discussed in the conclusion. 

We model the conduction electrons of our mesoscopic 
system in the absence of the core hole by the non- 
interacting Ho = a „Ci^a, whcre c\ ^ creates a 
particle with spin cr = ± in the orbital Lpi{v) {i ~ 
0, . . . ,N ~ 1). We furthermore assume that the perturb- 
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ing potential is a contact potential Vc = Vvc\rc){rc\, with 
Tc the location of the core hole and V the volume of the 
system. The diagonal form of the perturbed Hamiltonian 
isH = Ho+Vc = J2i,o- ^i4,a^i,<y^ where cJ,o- creates a par- 
ticle in the perturbed orbital tpi (r) . For the Fermi energy 
in the middle of the conduction band, Vc and the mean 
level spacing d are related through Sq — aictan(TrVc/d) [ij 
{So is negative since the core potential is attractive). 

A remarkable property of a rank one perturbation such 
as a contact potential is that all the quantities of inter- 
est for the x-ray edge problem can be either expressed in 
terms of the {e} and {A} or taken as independent random 
variables. For instance, ignoring for now the spin vari- 
able, the overlap between the many body ground states 
with M particles $o and \l/o of Hq and H is l] 
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For chaotic mesoscopic systems, we may assume that 
the unperturbed energy levels follow random matrix the- 
ory fluctuations and the wave function intensities are 
Porter-Thomas distributed (in particular at r^) Un- 
der these hypotheses, the joint probability distribution 
P(| e|, {A}) was derived by Aleiner and Matveev j3| (see 
|l3| for a generalization to the non-rank-one case): 
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with the constraint ei_i < A^ < e.^. Here /3 = 1 (circu- 
lar orthogonal ensemble, COE) or /3 = 2 (circular uni- 
tary ensemble, CUE) for the time reversal symmetric or 
asymmetric case, respectively. In the middle of the band 
Vi = Vc- In our case, it turns out to be necessary to in- 
clude boundary effects 0| , which modify the phase shifts 
away from the band center. This can be done simply by 
using a variable Vi given by 
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for i G [0, (-/V — l)/2], and the analogous form for 
i £ [{N — 1)/2,N — 1]. In addition, at perturbation 
strengths \vc\ ^ O.Sd, the lowest perturbed level Aq re- 
quires special treatment as its mean shift (Ao — eo) = 
— A^(i/[exp(|(i/uc|) — 1] becomes large. This gives rise to 
a second band in the x-ray spectra 1] that we do not 
discuss here. 

Using a Metropolis algorithm, we can therefore gener- 
ate ensembles of ({e},{A}) with the proper joint distri- 
bution, Eq. 0, and using Eq. ^ build the distribution 
of overlap for any value of the parameters M, N, Vc, 
and (3. To quantify the role of fluctuations, we define for 
reasons of comparison the special case of equidistant un- 
perturbed levels {e}, and will refer to it as the uniform 
or "bulkhke" case with overlap |Ahp (oc N~^f/'^ ). 



Results of such a simulation are displayed in Figs.^and 
12 The probability distributions of the overlap P(|Ap) 
shown in Fig. ^a) differ considerably even though the 
phase shift at the Fermi energy, 6f = arctan(7rz;i=M/rf), 
is similar. In fact, as expected, we find that both P(|Ap) 
and the bulklike result |Abp depend independently on A'^, 
A/, and 5f- However, performing a scaling by |Abp, see 
Fig.n^b), causes all curves with the same i5_f to collapse. 
Thus, the fluctuations of the overlap depend only on the 
value of the phase shift at the Fermi energy, 6f- 

The small |Ap part of P(|Ap) can be accurately repro- 
duced, as in the cases studied in Refs. by incor- 
porating the fluctuations of only the levels closest to the 
Fermi energy. For instance, we may multiply |Af,p by the 
fluctuating factor [1/(1 — Sf/tt)]{Xm+i — eM)/(AM+i — 
Xm)i which we call the "range-1" approximation. The 
resulting distribution (Fig. I^J is already similar to the 
exact one. We found an analytical result for this approx- 
imation for both the COE and CUE cases by treating 
two fluctuating levels and reducing the remaining ones 
to a single Gaussian variable 0|; it is seen to coincide 
with the range-1 result over the whole range of |Ap. 

We can similarly deflne "range-n" approximations 
which include the fluctuations of n levels on each side of 
the Fermi energy. The range-2 approximation matches 
the exact one for large perturbation strengths (Fig. I^J, 
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FIG. 1: (Color online) Probability distribution P(|Ap) of the 
overlap |Ap for different parameter sets {A^, M, —Vc/d}. (a) 
Although all Sf ~ — 7r/2, the P(|Ap) are visibly different as 
a result of different bulklike overlaps \Abf (increasing with 
the legend entries, COE case), (b) The curves coincide when 



scaled with lAi, 



The inset shows the CUE result with a 



noticeably different P(|A| ) distribution for small |A| origi- 
nating from the different wavefunction statistics. 
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FIG. 2: (Color online) Exact and approximate probability 
distributions of the overlap for (a) weak {vc/d — —0.25) and 
(b) strong {vc/d = -10) perturbation (COE, = 100, Af = 
50). The exact P(|Ap) is increasingly well reproduced by 
the range-1 and range-2 approximation, indicating that the 
fluctuations originate from levels around the Fermi energy. 
The dashed line denotes the range-1 analytical result. 



but does not improve significantly on the range-1 approx- 
imation for weak perturbations. We have checked, how- 
ever, that even in this case, it is sufficient to treat a few 
levels around Ep exactly to obtain the correct P(|Ap). 
This retrospectively justifies our use of a full random ma- 
trix to model our chaotic systems - the levels beyond the 
Thouless energy, where such a model does not apply [l5l | , 
do not affect the overlap distribution. 

Turning now to the absorption spectra, we limit our 
discussion to the K-edge (s-symmetry of the core electron 
wavefunction) and approach the mesoscopic x-ray edge 
problem using Fermi's golden rule 0, 

A{^) = 27Th-' Y.\{^f\D\<^^,)\H{Ef - i^S - c) , (5) 
/ 

where the sum is taken over all perturbed final states 
connected to the unperturbed groundstate <i>g = 
n,T=± n^o^ ^^Po\e operator D. (|c) = 

cl\0) is the empty band plus a core electron of spin 
a = +.) We are interested in processes involving the 
core hole; thus, the dipole operator can be written as 

D — const (^WjC^j -i-Cc + h.c^ . Since we assume the 

core electron wavefunction and the local part of the con- 
duction electron wavefunction are both s wave, Wj is 
related to the derivative of the perturbed orbital ipj in 
the direction e of the polarization of the x ray through 
Wj^e- VV'j(rc). 

First consider conduction electrons with the same spin 
as the excited core electron, and assume uj = wth so that 



the only possible final state is ^~jo = YijLo + ■ Without 
a perturbing potential, the only contribution is the direct 
process wmCm jf-Co- In the presence of a perturbation, 
however, the new and old orbitals are not identical, and 
terms with j < M (called replacement processes) also 
contribute coherently, giving 
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where is defined by generalizing Eq. ^ with level 

i (< M) replaced by M. Since for chaotic systems the 
derivative of the wavefunction, Vk~^ x jVg-V'jP, is known 
to have Porter-Thomas fluctuations uncorrelated with 
the wavefunction itself we can proceed as for the 
overlap to construct the distribution of K^I^yo |D|<i)g)p. 

Away from threshold, part of the x-ray energy can 
excite additional electrons above the Fermi energy in 
so-called shake-up processes. Their contribution is a 
straightforward generalization of Eq. I^. The number of 
these processes grows in principle exponentially with the 
energy of the x ray. However, it remains finite even for 
large lo — wth (and not very large for the values of M and 
N here) if one considers only shake-up processes that con- 
tribute significantly to absorption. Electrons with spin 
opposite to the excited core electron (spectator spin) af- 
fect the absorption only through AOC. 

Mesoscopic fiuctuations occur in both the photoab- 
sorption A{uj) and the excess energy uj — ujth] the for- 
mer arise from fluctuations in the dipole matrix elements 
and all energy levels, the latter from the levels above 
the Fermi energy. Here we limit ourselves to studying 
the average values of A{uj); fluctuations will be discussed 
elsewhere Q. We give {A{u!)) at the average excitation 
energy for a given final state, and measure all energies 
in mean level spacings d from the threshold Wth at which 
the core electron is excited to the perturbed level Am 
just above the Fermi energy. The normalized results are 
in Fig. |3|[a). The potential strength Vc = —10c? is cho- 
sen so that the Friedel sum rule is fulfilled {Sq « —tt/2, 
Si>o = 0). Shown are the spectra for the two spin types 
and the total result. We find a peaked edge for the opti- 
cally active spin sector as well as for the full spin result 
(obtained after convolution with the spectator spin) |0| . 

It is interesting to compare the chaotic case to the bulk- 
like situation (equidistant {e} and constant dipole matrix 
elements). Here, however, since each orbital channel acts 
"independently", one needs to specify whether the op- 
tically active channel is the one affected by the pertur- 
bation. As mentioned before, we have assumed s-type 
local symmetries, corresponding to a K-edge (^o — 1, 
<^/>o — 0, (5/=o ~ ~^/'^)- Since then only AOC plays a 
role, the bulklike situation gives a rounded edge; in con- 
trast, a peaked edge is found in the mesoscopic (chaotic) 
case [Fig. |3Jb)]. This striking difference is an effect of 
the coherent confinement in the chaotic system, in which 
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lie nanoparticles may be met within a few years. Another 
reahzation would be a double quantum dot with a tiny 
constriction Feasible with currently standard semi- 
conductor nanotechnology is an experiment based on a 
quantum dot array where the role of the core electron 
is taken by a localized electron bound to a (suitable) 
impurity level in the band gap. The excitation energy is 
provided by a micrometer laser that allows energy resolu- 
tion well below a mean level spacing. Our results predict 
a (slightly) peaked K-edge in the average photoabsorp- 
tion that becomes rounded for a bulk two-dimensional 
electron gas. 

We thank K. Matveev for several helpful discussions 
and I. Aleiner, P. Fuoss, Y. Gefen, I. Lerner, E. Muc- 
ciolo, U. Rofiler, I. Smolyarenko, W. Wegscheider, and 
D. Weiss for useful conversations. M.H. thanks the Hum- 
boldt Foundation for support. This work was supported 
in part by the NSF (DMR-0103003). 



FIG. 3: (Color online) X-ray photoabsorption spectra for 
(a) the mesoscopic (chaotic) and (b) the bulklike (non- 
fluctuating) case in terms of averaged cross sections (COE, 
iV = 100, M = 50,Vc/d = -10). (a) The mesoscopic spec- 
trum (flUed circles) shows a (slightly) peaked edge. Open 
circles indicate the individual spin contributions for the opti- 
cally active (full line) and spectator (dashed line) spins. The 
former arises as the sum of direct and replacement (dashed- 
dotted), as well as shake-up (dashed, crosses), processes. For 
comparison we give the bare spectrum (thin line) without ei- 
ther AOC or many-body effects, and the result with AOC 
only (direct process, thin dashed line), (b) Comparison of 
the peaked mesoscopic spectrum (circles) to the bulklike re- 
sults for the same optical channel (triangles, Iq — 1) that 
shows a rounded K-edge. The mesoscopic peak is, however, 
less pronounced than the bulklike L-edge (squares, Zo = 0). 

coupling is to the derivative of the wave function and so 
independent of the wave function itself. 

This peak is less pronounced, however, than that for 
the bulklike L-edge case (Iq — 0, so the optically excited 
and perturbed channels coincide). Because the dipole 
matrix elements Wi fluctuate in sign and magnitude in 
the mesoscopic case, the coherence of the replacement 
processes is less effective. Another difference between 
these two cases concerns their A^-dependence: whereas 
the peak sharpens with increasing A^ for the bulklike case, 
it diminishes with A^ in the chaotic situation. 

The basic physics we study - Fermi edge singulari- 
ties due to a localized (rank one) perturbation acting 
on a finite number of chaotic electrons via Fermi golden 
rule (dipole) matrix elements - is very general and al- 
lows a multitude of experimental realizations 0, 0, ^3 ■ 
The technical requirements for a direct implementation of 
mesoscopic x-ray absorption considered here using metal- 
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